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EXECUTIVE  SUMMARY 

A  program  of  basic  research  was  carried  out,  whose  objective  was  to  make  the  promising 
qualitative  thermal  wave  imaging  NDI  technique  a  truly  quantitative  tool,  which  can  have 
a  major  impact  on  the  rapid,  wide-area  inspection  of  Air  Force  aircraft  for  hidden 
corrosion.  Such  thermal  wave  images  of  a  square  foot  (or  so)  of  aircraft  skin  are  acquired 
on  short  time  scales  (a  few  seconds).  Techniques  were  developed  and  studied  to  measure 
skin  thickness  and  to  identify  corrosion  products.  The  resulting  corrosion  analysis  has 
been  shown  to  have  a  quantitative  capability  of  measuring  material  losses  as  small  as  1% 
on  aluminum  aircraft  skins,  and  can  be  implemented  rapidly  to  evaluate  regions  of 
corrosion  which  have  been  identified  in  the  thermal  wave  images.  The  technique  has  also 
been  shown  to  be  capable  of  measuring  intergranular  corrosion  in  the  vicinity  of  fasteners 
in  KC-135  and  B-52  wing  skins. 
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Introduction 

At  the  AFOSR  Logistics  Workshop:  Basic  Re  search  for  Logistics,  held  in  Dayton,  Ohio,  October  31- 
November  1, 1995,  Brig.  General  Michael  Moffitt  pointed  out  the  increased  reliance  of  the  Air  Force 
on  sustainment  of  aircraft,  as  the  result  of  a  decrease  in  the  purchases  of  new  aircraft.  He  noted  that 
corrosion,  fatigue,  and  cracking  pose  the  greatest  danger  to  the  sustainment,  pointing  out  that  22%  of 
the  C/KC-135  fleet  currently  is  in  depot  status,  as  compared  to  a  normal  expectation  of  10-15%.  One 
of  the  organizers  of  the  Workshop,  Dr.  Neal  Glassman  of  AFOSR,  next  reported  that  the  AFOSR 
Spring  Initiative  Review  had  identified  the  resulting  Logistics  issues  as  major  research  thrusts  for  new 
basic  research  grants.  The  research  conducted  under  AFOSR  Grant  No.  F49620-96-1-0166  was 
carried  out  in  response  to  some  research  issues  which  were  discussed  extensively  throughout  the 
ensuing  Workshop. 

At  the  outset  of  the  Grant,  thermal  wave  imaging  had  been  demonstrated  to  have  good  capability 
for  detection  of  hidden  corrosion,  as  well  as  of  disbonded  tear  straps  and  doublers.  However,  the 
quantitative  assessment  of  the  degree  of  corrosion  -  a  very  important  goal  for  the  successful  application 
of  the  technique  to  the  airworthiness  program  -  required  further  basic  research.  A  three-year  program 
of  research  had  been  proposed,  but,  because  of  budgetary  limitations,  the  research  was  funded  for 
eighteen  months  under  AFOSR  Grant  No.  F49620-96-1-0166.  This  Final  Report  presents  the  results 
of  that  eighteen  month  program  of  research. 
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Description  of  the  Experimental  Technique 

The  pulse-echo  thermal  wave  imaging  system,  which  was  developed  by  Wayne  State 
University1,  consists  of  a  pulsed  heat  source  (typically  high-power  photographic  flash  lamps),  an  IR 
video  camera,  and  image  processing  hardware  and  software,  all  of  which  is  controlled  by  a  personal 
computer. 

The  energy  from  the  pulsed  flashlamps  is  absorbed  at  the  surface  of  the  aircraft,  and  launches  a 
thermal  wave  pulse  into  its  skin.  When  this  pulse  is  reflected,  either  from  the  rear  surface  of  the  skin 
or  from  a  locally  corroded  interface,  the  reflected  portion  propagates  back  to  the  surface,  where  it 
modifies  the  time  dependence  of  the  temperature.  The  modification  is  greater  over  a  corroded  region 
than  over  an  uncorroded  one.  Time  delays  of  the  reflected  pulses  are  determined  by  the  thermal  wave 
transit  times  to  the  defect  and  back.  These  transit  times  to  the  rear  of  the  skin  and  back  are  shorter  for 
corroded  (thinner)  skin  than  for  the  uncorroded  skin.  The  returning  thermal  wave  reflections  are 
detected  by  means  of  the  IR  video  camera,  which  monitors  the  time-dependent  surface  distribution  of 
the  IR  emission  from  the  surface.  The  signal  from  the  camera  is  processed  in  real  time  by  fast  imaging 
hardware  and  software  in  the  computer,  and  the  thermal  wave  image  of  subsurface  corrosion  is 
displayed  on  a  video  monitor.  A  photograph  of  the  camera,  flashlamps  and  shroud  typical  of  the 
system  is  shown  in  Fig.  1 


Fig.  1  Photographs  of  the  exterior  (left)  and  interior  (right)  of  the  thermal  wave  imaging  head, 
showing  the  IR  focal  plane  array  camera,  aluminum  box,  and  two  linear  flash  lamps.  The  opening  of  the 
shroud  is  approximately  30  cm  on  a  side. 
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Description  of  the  Scientific  and  Technological  Barriers  for  Quantitative  Corrosion 
Thinning  Determinations,  Using  Thermal  Wave  Techniques 

Pulsed  thermal  waves  result  from  the  sudden  application  of  a  short  pulse  of  heat  to  the  surface 

of  a  solid  object.  This  launches  what  is,  in  effect,  a  temperature  wave  into  the  solid.  When  observed 

at  a  distance,  x,  beneath  a  planar  surface  of  the  object,  the  pulse  has  a  Gaussian  shape  given  by 


T(x,t)  = 


X 

4at 


(4ttat)1/2 

where  T  is  the  temperature,  t  is  the  time,  a  is  the  thermal  diffusivity  of  the  object,  and  A  is  a 
proportionality  constant  which  depends  on  the  heat  input.  At  a  defect  interface,  a  second  (reflected) 
pulse  is  generated  with  the  same  Gaussian  shape,  but  with  an  amplitude  which  depends  on  the 
reflection  coefficient  at  the  interface.  The  reflection  coefficient  is  a  function  of  the  ratios  of  the  thermal 
properties  on  the  two  sides  of  the  interface,  and  can  be  either  positive  or  negative.  The  reflected  pulse 
appears  to  have  originated  at  the  image  of  the  object  surface,  as  seen  in  the  "mirror"  of  the  defect 
interface.  Thus,  it  looks  like  a  thermal  wave  "echo"  from  that  interface.  When  it  reaches  the  surface 
(x=0),  it  modifies  the  time-dependence  of  the  surface  temperature  by  contributing  an  additional  term  to 
the  inverse  square  root  of  time  behavior  of  Eq.  (1)  at  x=0: 


(1) 


T(0,t) 


( 


A 

(4mt)1/2 


1  +  Re 


(2d)2  ^ 
4at 


(2) 


V  J 

Here,  d  is  the  depth  of  the  defect,  and  R  is  the  reflection  coefficient.  The  time  at  which  the  additional 
term  in  Eq.  (2)  becomes  significant  depends  on  d2/a.  Thus,  one  can,  in  principle,  measure  the  depth 
(e.g.,  thickness  of  a  corroded  skin)  by  observing  the  time  at  which  the  deviation  from  Eq.  (1)  occurs. 
However,  this  measurement  is  complicated  by  variations  in  the  (unknown)  value  of  R,  and  by  the 
effects  of  three-dimensional  diffusion,  which  have  been  ignored  in  this  simplified  one-dimensional 
picture.  Also,  if  one  waits  a  bit  longer,  one  gets  additional  terms  in  Eq.  (2),  due  to  multiple  reflections 
of  the  pulse.  It  should  also  be  kept  in  mind  that  thermal  waves  are  extremely  dispersive,  so  that  the 
pulses  broaden  dramatically  as  time  progresses,  and  the  multiple  "echoes"  all  overlap.  Nevertheless,  if 
the  surface  temperature  is  observed  at  short  enough  times  after  the  flash,  one  can  essentially  eliminate 
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the  effects  of  both  multiple  reflections  and  three-dimensional  diffusion,  and  make  a  determination  of 
the  "time-of-flight".  A  mathematical  model2  based  on  three-dimensional  wave  scattering  theory, 
accurately  predicts  the  temperature-time  profiles  at  a  given  point,  and  also  predicts  the  variations  of  the 
temperature  on  the  surface  above  an  arbitrarily-shaped  defect.  This  theory  has  been  used  to  calculate 
theoretical  images  of  defects  with  several  simple  geometries  in  steel  and  composite  materials,  and  has 
been  confirmed  by  comparison  with  experimental  thermal  wave  images  of  corresponding  test 
specimens. 

The  challenge  for  this  research  project,  both  scientifically  and  experimentally,  was  to  seek  ways 
to  exploit  the  early-time  behavior  so  as  to  determine  the  nature  of  the  corrosion  and  to  measure  the 
residual  skin  thickness.  Since  the  characteristic  transit  times  for  aluminum  skins  are  very  short 
(comparable  with  the  frame  times  of  the  IR  imaging  systems  which  the  were  state-of-the-art  at  the 
outset  of  the  project),  we  considered  alternative  strategies  for  processing  the  data  from  such  cameras. 
We  also  considered  the  potential  benefits  of  utilizing  more  advanced  IR  focal  plane  arrays  which 
became  available  during  the  period  of  performance  of  the  Grant.  Both  approaches  proved  to  be 
beneficial. 

Data  Processing:  Implementation  of  a  Corrosion  Analysis  Program  in  the  Thermal 
Wave  Imaging  Software 

Another  practical  consideration  which  was  noted  in  our  Research  Proposal  on  this  Grant  is  that 
an  aircraft  skin  provides  no  thick  background  area  for  determining  image  contrast.  Test  specimens 
which  had  been  used  in  our  preliminary  research  had  been  fabricated  from  thick  metal,  which  made  a 
convenient  background  reference  for  pulse-echo  analysis.  However,  in  an  aircraft  skin,  the  reference 
will  have  to  be  the  undamaged  skin  areas,  which  are  of  the  order  of  1mm  in  thickness.  This  reference, 
must  then  be  compared  with  the  corroded  areas,  which  may  be  only  a  few  percent  thinner. 
Nevertheless,  despite  the  lack  of  a  thick  reference,  we  noted  in  our  Proposal  that  we  had  achieved 
sensitivity  to  corrosion  thinning  of  a  few  percent  material  loss.  The  problem  was  to  determine  a 
strategy  for  making  quantitative  measurements  of  this  percentage  loss.  We  developed  such  a  strategy, 
based  upon  the  following  analysis,  described  in  greater  detail  in  the  Ph.D.  Dissertation  of  Dr.  Xiaoyan 


9 


Thermal  Wave  Imaging  of  Hidden  Corrosion  FINAL  REPORT 

R.L.  Thomas,  L.D.  Favro  and  P.K.  Kuo  WAYNE  STATE  UNIVERSITY 

AFOSR  Grant  No.  F49620-96-1-0166 


Han  of  Wayne  State  University.3 


Consider  a  sample  which  has  two  insulated  surfaces  Si,  S2  at  z  =  0  and  z  =  d  shown  in  Fig.  2. 
For  one-dimensional  heat  propagation  along  the  z-direction  with  the  initial  condition, 


T(z,0+)  =  8(z)  ,  (3) 

Applying  Eq.  (2)  to  the  case  of  a  thermally  thin  material,  and  setting  R=1  (insulating  boundaries),  we 
have 


(z-2nd)2 


T(z,t)  = 


4at 


(47tat) 


1/2 


(4) 


Fig.  2  Schematic  diagram  of  a  perfect  sample  with  two  insulated  surfaces  Si,  S2  at 
z  =  0  and  z  =  d  respectively. 

This  solution  is  good  for  theoretical  calculations  when  we  deal  with  short  time  scales  since  it  converges 
quickly  for  small  t.  However,  when  we  need  to  look  at  longer  time  scales,  we  need  another  form  of  the 
solution  to  the  heat  equation, 

92T(z,t)  1  8T(z,t) 

d2z  a  3t  ' 

Choosing  the  solution  in  the  form, 


T(z,t)  ~  cos(kz)e  71  , 


(6) 
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after  some  manipulation3,  we  obtain  an  expression  of  the  form 


oo  i*u 

T(z,t)  =  ^Ancos(^z)e  d: 

£o  d 


(7) 


According  to  the  boundary  condition  expressed  by  Eq.  3,  we  have. 


jr  Ancos(^z)  =  5(z)  . 

n=0  d 

(8) 

± a4>cos(Hz)cos(Mz)  =  Jdz8(z)cos(^pz)  , 

(9) 

^  AnJ  dzcos(— z)  =  Jdz5(z)  ,  and 

n=0  0  d  0 

(10) 

Eqs.  8-10  lead  to 

2 

Am  =  -  ,  for  m^O  , 
d 

(11) 

and. 

|  -a 

II 

o 

< 

(12) 

So,  the  general  solution  to  Eq.  5  can  be  written  as. 

2  2 

1  00  nK  an2n2 1 

T(z,t)  =  -[l  +  2£cos(— -  z)e  <*2  ]  . 
d  d 

(13) 

At  the  surface  z  =  0,  the  temperature  can  be  expressed  as, 

,  „  anV  ( 

T(0,t)  =  -[l  +  2^e  d*  ]  . 

d  n=l 

(14) 

«o  an2*2 1 

Finally,  at  a  time  t,  when  2^T  e  d2  «  1  , 

n=l 

(15) 

we  have, 
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T(0,t)ocl  .  (16) 

d 

At  sufficiently  long  time,  t,  Eq.  16  can  be  used  to  measure  the  skin  thickness  d.  However,  if  t  is  too 
long,  the  assumed  boundary  conditions  will  be  poorer  approximations,  and  will  lead  to  systematic 
errors.  Thus,  it  is  important  to  wait  long  enough  after  the  flash  so  that  the  sample  achieves  local  thermal 
equilibrium,  but  not  so  long  that  heat  losses  to  the  environment  or  to  thicker  regions  of  the  sample 
become  important.  As  a  practical  matter,  for  typical  aluminum  aircraft  skins,  the  term  represented  by 
Eq.  15  is  negligibly  small  after  about  20ms  following  the  flash.  In  our  experiments,  we  typically 
acquire  a  first  image  at  about  80ms  following  the  flash,  so  we  expect  the  simple  approximation  of  Eq.  16 
to  hold,  and  utilize  this  approximation  in  an  algorithm  to  determine  the  fractional  loss  in  thickness  by  the 
expression 

AT2  ~  AT,  _  dt  —  d2 
AT2  d, 

We  have  implemented  such  an  algorithm,  based  on  Eq.  (17),  in  our  thermal  wave  imaging  corrosion 
analysis  window  software,  and  have  tested  it  on  a  number  of  intentionally  corroded  regions  of  different 
thicknesses  of  aluminum  skins.4  The  window,  shown  in  Fig.  3  (left),  contains  “draggable”  regions  of 
interest  for  corrosion  analysis,  including  some  which  are  placed  over  uncorroded  reference  regions  of 
the  specimen.  The  percentage  corrosion  of  the  other  regions  of  interest  read  out  immediately  as  they  are 
dragged  above  the  region  of  the  image  containing  suspected  rear-surface  corrosion.  The  results  of  the 
read-out  of  percentage  corrosion  are  plotted  on  the  right  in  Fig.  3,  and  indicate  excellent  agreement  with 
direct  micrometer  measurements  of  percentage  material  loss  for  greater  than  1%  corrosion  material  loss. 
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Fig.  3  Corrosion  analysis  window  (left),  containing  a  thermal  wave  image  of  an  aluminum  test  panel, 
the  rear  surface  of  which  has  five  intentionally  corroded  regions.  Also  shown  in  the  image  are  smaller 
squares  which  indicate  the  region  being  analyzed  and  two  reference  regions.  The  resulting  comparison 
with  micrometer  measurements  shown  on  the  right  for  this  and  similar  panels,  indicates  excellent 
agreement  from  less  than  one  percent  to  nearly  thirty  percent  material  loss. 


Using  the  same  algorithm,  we  can  also  “paint”  regions  of  corrosion  identified  on  a  thermal  wave 
image,  “painting”  regions  whose  material  loss  exceeds  a  threshold  value  which  can  be  chosen  by  the 


Loss  >  15%  ;  Area  1.17  sq.in  Loss  >  25%  ;  Area  0.28  sq.in 


Fig.  4  (Left)  a  raw  thermal  wave  image  of  a  corroded  riveted  two-layer  sample  which  contains  rear 
surface  corrosion  of  the  front  layer  with  the  corrosion  products  in  place,  supplied  by  Boeing. 
(Right)  four  analysis  images,  prepared  using  our  real-time  algorithms.  In  red  are  shown  the 
areas  with:  >5%  (top  left),  >10%  (top  right),  >15%  (lower  left),  and  >20%  (lower  right) 
corrosion  material  loss.  Also  printed  out  as  part  of  the  analysis  are  the  actual  areas  (in  square 
inches)  of  corrosion  exceeding  these  bounds. 
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Advanced  IR  Focal  Plane  Array  Imaging  System 

A  camera,  whose  technology  was  based  on  our  one-of-a-kind,  512x512  Imager,  which  was 
described  in  detail  in  the  Proposal,  since  became  available  as  a  commercial  product  and  with  a  price 
that  is  competitive  with  other  FPA  cameras.  This  camera  (Amber/Raytheon  Galileo)  is  a  256x256 
array,  with  a  base  frame  rate  of  125  Hz,  and,  like  our  512x512  camera,  operates  in  the  snapshot  mode. 
This  camera  also  has  a  modified  version  of  the  windowing  capability  which  we  had  designed  into  our 
research  camera,  and  reaches  a  maximum  frame  rate  of  2kHz.  In  addition,  whereas  the  512x512 
research  camera  has  liquid  nitrogen  cooling,  this  commercial  camera  has  its  own  Stirling  Engine 
cooler,  making  it  quite  suitable  for  use  in  a  hangar  environment.  In  conjunction  with  our  technology 
transfer  effort,  and  with  funds  from  other  (non-AFOSR)  sources,  we  have  purchased  one  of  these 
commercial  cameras,  and  utilized  it  in  support  of  this  research.  Because  of  the  fact  that  the  Research 
Program  terminated  at  eighteen  months,  rather  than  thirty  six  months  effort,  the  Galileo  camera  was 
utilized  less  extensively  than  originally  planned.  However,  we  initiated  its  use  to  obtain 
temperature/time  curves  at  the  higher  frame  rates,  and  anticipate  that  this  approach  could  prove  to  be  a 
useful  complement  to  the  corrosion  analysis  algorithm  described  above. 

Technology  Transfer  Project:  KC-135  and  B-52  Wingfastener  Corrosion 

At  the  request  of  Oklahoma  City  ALC,  ARINC  conducted  a  round-robin  study  of  the  capability 
of  various  NDE  techniques  to  detect  intergranular  corrosion  in  the  vicinity  of  steel  fasteners  in  KC-135 
and  B-52  wingskins.  We  participated  in  the  study  January  27, 1997,  which  included  a  blind  test  of  80 
such  fasteners.  Example  images  of  four  of  these  fasteners  are  shown  in  Fig.  5.  We  had  zero  “false 
calls”  in  our  blind  test,  and  correctly  identified  87.5%  of  the  corroded  fasteners  at  a  rate  of  0.5 
fasteners  inspected  per  minute.  Thermal  wave  imaging  is  a  very  promising  inspection  technique  for 
this  important  ALC  inspection  requirement. 
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Fig.  5  Images  of  an  uncorroded  (top  left)  and  three  corroded  wing  fastener  regions,  taken  from  an 
ARINC  test  of  corroded  C/KC-135  wing  skins.  The  intergranular  corrosion  shows  up  as  an  irregular 
boundary  around  the  image  of  the  fastener. 
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